Hubbard model versus t-J model: The one-particle spectrum 
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The origin of the apparent discrepancies between the one-particle spectra of the Hubbard and t-J 
models is revealed: Wavefunction corrections, in addition to the three-site terms, should supplement 
the bare t-J model. In this way a quantitative agreement between the two models is obtained, 
even for the intermediate- (7 values appropriate for the high-T c cuprate superconductors. Numerical 
results for clusters of up to 20 sites are presented. The momentum dependence of the observed 
intensities in the photoemission spectra of S^CuC^Ch are well described by this complete strong- 
coupling approach. 
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The past ten years the interest in Hamiltonians mod- 
eling correlated electrons has greatly increased. This is 
largely due to the mysterious normal-state properties of 
the high-T c cuprate superconductors and the realization 
that local correlations on the Cu sites are large in these 
compounds. The two most intensively studied Hamilto- 
nians are the Hubbard and t-J models. Although at first 
sight such models may seem oversimplified, it has been 
argued that the two dimensional square-lattice Hubbard 
model captures the essential dynamics of the cuprates 
in a window of a few eV around the chemical poten- 
tial The materials belong to the intermediate-// 
regime, U/t » 10, t » 0.3 - 0.4 eV [§,[§. Further support 
for a single-band Hubbard approach comes from mea- 
sured one-particle and optical spectra of the cuprates. 
These show a dramatic reshuffling of intensity as a func- 
tion of doping on the scale of a few eV |^|| . This char- 
acteristic behavior is a fingerprint of strong correlations 
and is described well by the Hubbard model • 

Although the t-J model is frequently referred to as be- 
ing the large-?/ limit of the single-band Hubbard model, 
the (one-particle) spectra of the two models differ sig- 
nificantly in the intermediate-?/ regime of interest. (For 
a recent review on numerical results, see Ref. |Mj) Prob- 
ably because of this, the two models are often treated 
as being "distinct" . One obvious difference concerns 
the momentum occupation number at half filling, 
equal to the one-electron removal sum rule. For the t-J 
model this quantity is trivially equal to 1/2, while for the 
intermediate-?/ Hubbard model it is a strongly varying 
function of k. Even though in the infinite- U limit 
tends to 1/2, as in the t-J model, the convergence to this 
value by increasing U is very slow. 

In this Letter we show that these discrepancies can be 
attributed to the neglect of certain terms in the large- 
U perturbation theory. Specifically, the weight distribu- 
tions in the spectral functions show a high sensitivity to 



the first-order renormalization of the wave functions. The 
complete strong-coupling approach is found to quantita- 
tively reproduce the Hubbard spectra for U > W, where 
W is the U — band width. For these U values the per- 
turbation series in t/U agrees well with the full Hubbard 
Hamiltonian. 

The derivation of the strong-coupling Hamiltonian by 
means of a canonical transformation, written as an ex- 
pansion in t/U, is well known (^^]. Starting point is the 
Hubbard model, 



H = V 
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Here 71,-,-f is the occupation number of the original Hub- 
bard fermions dj )CT , and 8 is a vector connecting nearest 
neighbors. After the transformation the states in the 
lowest Hubbard sector (no doubly occupied sites) are de- 
scribed to order t 2 /U by the strong-coupling Hamilto- 
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where J = 4t 2 /U. The motion of the above transformed 
c,', CT fermions will leave the number of doubly occupied 
sites unchanged, and Si and ni jCT are the spin and occu- 
pation number of these new fermions. Simply neglecting 
the 'three-site term', i.e. the second line of (||), one ob- 
tains the t-J model. 

Why the full strong-coupling model has received far 
less attention than the t-J model is unclear. It might be 
expected that the three-site term is very important in the 
one-hole limit (and low doping region) since, unlike the 
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nearest-neighbor hopping, a direct next-nearest-neighbor 
hopping process does not frustrate the spin system. How- 
ever, it was shown by a numerical study ]l0[ , spin-polaron 
approaches ]ll]Jl2]] and a linear-spin-wave self-consistent 
Born (LSW-SCB) calculation Jl3| that the main effect of 
this term on the quasiparticle dispersion is an increase of 
the bandwidth of the dominant 'quasiparticle peak'. 

By virtue of its construction via a canonical transfor- 
mation, the Hamiltonian (Q) reproduces the eigenvalue 
spectrum of the Hubbard model for large U. To obtain 
further information, such as matrix elements or expec- 
tation values, operators acting in the Hilbert space of 
the Hubbard model have to be subjected to the same 
canonical transformation. For the electron annihilation 
operator this procedure leads to the expression Mm , 



Pa t .„P = Pci,„P 



p. 



(3) 



P is defined as the projector on states with no dou- 
ble occupancy of the transformed fermions (states with 
(V) = 0). The first-order correction terms (second line) 
describe the dressing of the annihilation operator with 
'virtual' charge fluctuations involving doubly-occupied 
sites. While their prefactor t/U seems to suggest that 
they represent only a minor correction, their impact on 
the calculated spectra is dramatic - one reason being 
that the summation over neighbors introduces the co- 
ordination number z as an additional 'hidden' prefac- 
tor. When comparing strong-coupling calculations di- 
rectly with experiments it is therefore crucial to take 
these corrections into account. For the optical conduc- 
tivity such corrections were shown to lead to qualita- 
tively different results Another consequence is a 
rapid reshuffling of weight between the upper (UHB) and 
lower-Hubbard band (LHB) as a function of hole doping, 
observed in both one-particle and optical experiments. 
The first-order corrections to transformed operators such 
as Eq. (||) lead to modified expressions for sum rules that 
consistently describe the rapid intensity changes observed 
in experiments ||||,|]]9) ■ 

For the numerical calculations of the one-particle spec- 
tra we implemented the strong-coupling model Eq. (2). 
Electron removal is described by the operator Eq. (3). 
This operator does not change double-occupancy and 
only the LHB part of the Hubbard spectrum is pro- 
duced. (There is, in fact, a similar independent proce- 
dure to calculate the UHB response as well.) The t/U 
term in Eq. (j^) amounts to the removal of an electron 
from a neighbor of i. This leads to an extra phase factor 
exp(zk-<5). Because of this the total weight is enhanced at 
r and reduced at (ir, ir) and Eq. (^) leads to the correct 
behavior of n(k) (see Fig. 2 of Ref. [§). Note that the 
implementation of the first-order correction is straight- 
forward and does not increase memory or CPU time in 
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FIG. 1. The one-particle spectrum of the two-dimensional 
4x4 strong-coupling model (solid line) compared with the 
Hubbard model (dotted line), at half filling (left) and for two 
holes (right). U/t = 10. Only the LHB part of the Hubbard 
spectrum is shown. Dashed line is the chemical potential in 
the two-hole case. Hubbard data reproduced from Ref. 

a Lanczos-like diagonalization calculation. In a similar 
way such operator corrections can be taken into account 
in analytic approaches, for instance in the LSW-SCB cal- 
culation of the one-hole spectrum. 

In Fig. |l| we compare the two-dimensional 16-site 4x4 
strong-coupling model with the LHB electron removal 
and addition spectra of the Hubbard model, reproduced 
from Ref. {nj. Note that U — 10 is comparable to the 
bandwidth W = 8. No shift or scaling has been applied 
to the spectra, showing that the value of the gap between 
the UHB (not shown, starting at energy w 8i) and LHB 
is accurately described. (While there is no UHB in the 
strong coupling model, the magnitude of the energy gap 
still can be inferred from the relation A = U — 2 • E m i n , 
with E m i n w 2t the minimal ionization energy of the 
model.) The resulting gap of « 6t is consistent with the 
measured insulating gap of « 2eV. 

As can be seen, excitation energies and weights of the 
dominant features in the spectra are reproduced correctly 
by the strong-coupling approach. The agreement is best 
for the physically most relevant low excitation energies, 
but the strong coupling approach also gives a complete 
description of the higher excitations, with energies not 
much smaller than U. The total width is slightly overes- 
timated by the second-order t 2 /U terms in the Hamilto- 
nian, due especially to those peaks which obtain a large 
energy correction by the three-site terms, for instance at 
(•7T, 7r). There is no significant deterioration of the agree- 
ment upon doping, demonstrating that the discrepancies 
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between the doping dependences of low energy sum-rules 
for t — J and Hubbard originate from matrix element ef- 
fects, rather than an incompleteness of the Hilbert space 
of the strong coupling model. We explicitly checked that 
second-order corrections to the Fermi operators lead to 
only minor changes of the peak intensities. For U values 
below 10 the agreement between Hubbard and its strong- 
coupling limit gradually deteriorates, which we interpret 
as a breaking down of the perturbation approach. 

Stimulated by the good agreement in the 4x4 system 
we now present large-U Hubbard spectra for clusters up 
to 20 sites. Because the momenta of the 16 and 18 site 
clusters nicely complement each other and are all on the 
high-symmetry axis, we show a mix of these two clus- 
ters in Figs. 2 and 3. The curves on the left show a scan 
from (0, 0) to (tt, tt). The plots on the right go from (0, 0) 
to (tt, 0) and from (tt, 0) to (tt,tt). We shifted all 4 x 4 
spectra by the same small energy such that the lowest 
peak at k = (0, 0) is at the same energy for both clus- 
ters. For the two-hole case we shifted the t-J data by an 
additional 0.3i to have the chemical potentials at roughly 
the same position. The ground-state momentum is (0, 0) 
in all cases (there is an extra degeneracy for 16-sites, 2 
holes) . In both figures the results are compared with the 
"bare" t-J model. 

As mentioned before there are substantial differ- 
ences between t-J and strong coupling. At half-filling 
the low energy 'quasiparticle peaks' are at higher 
energies (implying a larger Hubbard gap) for the 
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FIG. 2. The momentum dependence of the elec- 

tron-removal spectrum of the strong-coupling Hubbard model 
(solid line) compared with the t-J model (dotted line). The 
plot is a mix of data from the 16 and 18 site clusters. 
U/t = 10. n = 1. The spectra for k = (tt/2,tt/2), (tt/2,0), 
(tt, 0) and (7r/2, tt) are calculated using the 2D 16-site cluster; 
the others, using the 18-site cluster. 
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FIG. 3. Combined electron addition and removal spectra 
for the two-hole ground state. See caption of Fig. ^. 

t-J model and their dispersion has a smaller width. 
Moreover, the momentum dependence of the weights is 
qualitatively different. The t-J spectrum at (0,0), and 
indeed through much of the Brillouin zone, is more com- 
pact than for large-£7 Hubbard. The strongest difference 
is at (tt, tt). Here strong coupling shows a little peak at 
high energy, while the t-J model has pronounced struc- 
tures at a more intermediate energy. This reflects the 
small value of for k = (jr, tt) in the Hubbard model, 
which is reproduced by the first-order corrections to the 
Fermi operator Eq. (||). 

The spectra at half filling can be compared with re- 
cent angle-resolved photoemission data on insulating 
Sr2Cu02Cl2 Jl(|. A significant feature of these measure- 
ments is the strong momentum dependence of the weights 
of the spectra. The region in momentum space with the 
most intense peaks is not centered around (7r/2,7r/2), 
but displaced towards (0,0). Beyond (tt/2,tt/2) the in- 
tensity drops very sharply and seems to vanish around 
(tt, tt) p~7| . As shown in Fig. 4 these features are well 
reproduced by the strong-coupling Hubbard model. In 
contrast, the intensities in t-J are much more symmetric 
around (tt/2,tt/2) and change only little in that part 
of momentum space. On average the peaks in strong- 
coupling Hubbard are more intense than those in the 
t-J model, at the expense of the incoherent part of the 
spectrum. Note that also in the t-J model there are 
differences in weight between k and k + (tt, tt). This is 
due to the coupling of the created hole to quantum-spin 
fluctuations in the ground state of the Heisenberg anti- 
ferromagnet Experimentally the strong changes at 
the antiferromagnetic Brillouin-zone boundary look like 
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FIG. 4. Weights (proportional to the radius of the circles) 
of the lowest energy peak for the strong-coupling Hubbard 
model (left) and t-J model (right). Data collected from the 
16, 18 and 20 site clusters. U/t = 10, n = 1. 

a kind of 'pseudo Fermi surface': The spectral weight 
observed drops significantly, but an actual Fermi-level 
crossing does not take place. We note that the large 
difference observed between (tt, 0) and (tt/2, 7r/2) in the 
photoemission measurements actually is in disagreement 
with the near degeneracy of these momenta in both the 
standard Hubbard and t-J model. As discussed by vari- 
ous authors this may be explained by adding extra next- 
neighbor hopping terms to the Hamiltonian fl"3|]l8| . 

The increased quasiparticle weight and the strong drop 
of intensity when passing the half-filled Fermi surface 
are consistent with a continuous cross-over to a small- U 
spin-density- wave (SDW) like shadow-band scenario [^9| . 
Surely Fig. ^ shows that such a strong drop in weight is 
not necessarily a weak-coupling feature. 

The interpretation of the doping dependence of the 
spectral function is a subtle issue pc|-|2^] , and will not be 
discussed here. For roughly 10% holes (Fig. ||) there are 
intense features appearing in the electron addition spec- 
trum for all momenta outside the magnetic zone. These 
features are more intense than for t-J and because of 
the sum rule the extra weight must come from the up- 
per Hubbard band. This weight-transfer effect has been 
clearly observed in O-ls X-ray absorption S. 

In conclusion, we have shown that the strong-coupling 
approach is a good alternative to Hubbard calculations 
even in the intermediate U regime relevant for the high- 
T c superconductors. Only the complete strong-coupling 
model gives a faithful representation of the properties 
of the Hubbard model. The corrections to t-J lead to 
a smooth crossover to small U. For comparisons with 
experiments such as angle-resolved photoemission, O-ls 
X-ray absorption (and optical conductivity) it is impor- 
tant to use the renormalized wavefunctions (rotated op- 
erators), as these can lead to even qualitatively different 
results. 
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